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Abstract 

Generic classically integrable boundary conditions for the affine Toda field the- 
ories (ATFT) are investigated. The present analysis rests primarily on the underlying 
algebra, defined by the classical version of the reflection equation. We use as a proto- 
type example the first non-trivial model of the hierarchy i.e. the ATFT, however 
our results may be generalized for any An^ (n > 1). We assume here two distinct 
types of boundary conditions called some times soliton preserving (SP), and soliton 
non-preserving (SNP) associated to two distinct algebras, i.e. the reflection algebra 
and the (q) twisted Yangian respectively. The boundary local integrals of motion are 
then systematically extracted from the asymptotic expansion of the associated transfer 
matrix. In the case of SNP boundary conditions we recover previously known results. 
The other type of boundary conditions (SP), associated to the reflection algebra, are 
novel in this context and lead to a different set of conserved quantities that depend 
on free boundary parameters. It also turns out that the number of local integrals of 
motions for SP boundary conditions is 'double' compared to those of the SNP case. 
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1 Introduction 



Integrability in the bulk has admittedly attracted a great deal of research interest in recent 
years, however after the seminal works of [H [21 [3] particular emphasis has been given on 
the issue of incorporating consistent boundary conditions in integrable models. This shed 
new light into the bulk theories themselves, and also opened the path to new mathematical 
concepts and physical applications. In a more general setting the investigation of both classi- 
cal and quantum integrable systems, particularly those with non-trivial boundary conditions, 
turns out to be quite significant especially after the recent advances within the AdS / CFT cor- 
respondence [I] uncovering the important role of integrability [5]. A crucial question within 
this frame is what would the physical implications be in both gauge and string theories once 
non-trivial consistent boundary conditions, especially the ones that may modify the bulk 
behavior, are imposed to the associated lattice and continuum integrable models (for some 
recent results see [6] and references therein). Therefore studies concerning the existence of 
consistent boundary conditions that preserve integrability are of particular significance and 
timeliness not only for the integrable systems themselves, but for other active research fields. 

The central purpose of the present article is the investigation of classical integrable models 
when general boundaries that preserve integrability are implemented. Among the various 
classes of integrable models we choose to consider here a particular class that is the afline 
Toda field theories (ATFT) [3, [8]. The prototype model of this class is the sine-Gordon model, 
which has been extensively studied both in the bulk [9] as well as in the presence of non- 
trivial integrable boundary conditions [TO]. Generic affine Toda field theories with classical 
integrable boundary conditions were first analyzed more than a decade ago in [IT] . A different 
point of view, although regarding the same class of boundary conditions^ analyzed in [11] , is 
presented in [12]. Specifically, in [12] the A^p ATFT with 'dynamical' boundary conditions 
-that is a quantum mechanical system is attached at the boundary — is investigated. Further 
studies regarding the boundary ATFT at both classical and quantum level may be also found 
in various articles (see e.g. [T3] [T8]). 

Although the analysis in [11] seems quite exhaustive it turns out that in simply-laced 
ATFT a whole class of consistent boundary conditions is absent. Our main objective here 
is to systematically search for all possible boundary conditions in An ATFT and eventually 
implement the missing ones. More precisely, we assume two distinct types of boundary 
conditions called soliton preserving (SP), and soliton non-preserving (SNP) associated to two 
distinct algebras, i.e. the reflection algebra [2] and the twisted Yangian [191 120] respectively 

2 by 'same class of boundary conditions' we mean that in both studies [111112] a common underlying algebra 
-(classical) q-twistcd Yangian- is implicitly assumed. Note however that the analysis in [IT] is classical while 
in [12] is quantum. 
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(see also relevant studies in [T8| I2T1 [22| 1231 12^ 125]). 

Depending on the choice of boundary conditions certain physical behavior is entailed. 
Specifically, in the context of imaginary A n ATFT the boundary conditions introduced in 
[TT] . known as SNP, oblige a soliton to reflect to an anti-soliton. In real A„ ATFT on the 
other hand such boundary conditions lead to the reflection of a fundamental particle to itself. 
Recall that fundamental particles in real ATFT are equivalent to the lightest bound states 
(breathers) of the imaginary theory provided that (3 —>■ i/3 ((3 is the coupling constant of 
the theory). It is however clear that another possibility arises, that is the implementation 
of certain boundary conditions that lead to the reflection of a soliton to itself in imaginary 
ATFT or to the reflection of a fundamental particle to its conjugate in real ATFT. These 
boundary conditions are known as soliton preserving and have been extensively analyzed in 
the frame of integrable quantum spin chains [21], [22], [26] |30j. 

Notwithstanding SP boundary conditions are somehow the obvious ones in the framework 
of integrable lattice models they have remained elusive in the context of An ^ ATFT for quite 
a long time. Note however that in quantum spin chains in addition to the well studied SP 
boundaries SNP boundary conditions were first introduced in [31] and further analyzed and 
generalized in [211 E2J EU [25]. It is thus our primary objective here to complete the study of 
integrable boundary conditions in ATFT by introducing and fully analyzing the novel (SP) 
boundary conditions. 

The outline of this article is as follows: in the next section we present the basic prelimi- 
nary notions regarding the algebraic setting for classical models on the full line and on the 
interval. In our analysis we adopt the line of attack described in e.g. [32], and in [331 El] 
for boundary systems. More precisely we introduce the classical Yang-Baxter equation and 
the underlying algebra for the system on the full line. In the situation of a system on the 
interval we distinguish two types of boundary conditions based on the classical versions of 
the reflection algebra (SP) and (q) twisted Yangian (SNP). Next the A n ATFT on the full 
line is reviewed and an explicit derivation of the local integrals of motion by solving the 
auxiliary linear problem [32] is presented. In section 3 being guided by the same logic and 
adopting Sklyanin's formulation [2] we rederive the integrals of motion of the A n ^ ATFT 
with SNP boundary conditions. Note that analogous strategy was followed in [33] and [31] 
for the classical boundary sine-Gordon and vector NLS models respectively. Our results are 
in agreement with the ones deduced in [TT]. In section 4 we introduce for the first time the 
novel boundary conditions (SP) within the context of ATFT. Explicit expressions of the as- 
sociated local integrals of motion are deduced from the asymptotic expansion of the classical 
transfer matrix. It is worth stressing that the induced integrals of motion depend on free 
boundary parameters as opposed to the SNP case. In the last section a discussion on the 
entailed results is presented and several directions for future investigations are proposed. 
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2 Preliminaries 



The analysis of the ATFT with integrable boundary conditions will rely on the solution of the 
so called auxiliary linear problem [32]. Before we proceed to the study of classical integrable 
models with consistent boundary conditions it will be instructive to recall the basic notions 
in the periodic case. Let $ be a solution of the following set of equations 



(2.1) 
(2.2) 



with U, V being in general nxn matrices with entries functions of complex valued fields, their 
derivatives, and the spectral parameter A. Compatibility conditions of the two differential 
equation (12. ip . (12. 2p lead to the zero curvature condition [351 



U- V' + 



U, V 



0. 



(2.3) 



The latter equations give rise to the corresponding classical equations of motion of the system 
under consideration. The monodromy matrix from (12. ip may be written as: 



T(x, y, A) = Vexp{ j U(x', t, \)dx' 



(2.4) 



with T(x,x, A) = 1. The monodromy matrix satisfies apparently (12.11) . and this will be 
extensively used in the present analysis. On the other hand within the Hamiltonian formalism 
the existence of the classical r-matrix, satisfying the classical Yang-Baxter equation 



ri 2 (Ai - A 2 ), ri 3 (Ai) + r 23 (A 2 ) + r 13 (Ai), r 23 (A 2 ) 



0, 



(2.5) 



guarantees the integrability of the classical system. Indeed, consider the operator T(x, y, A) 
satisfying 



Ti(x,y,t,\i), T 2 (x,y,t,X 2 ] 



r 12 (\i - A 2 ), T 1 (x,y,t,X 1 )T 2 (x,y,t,X 2 ) 



Making use of the latter equation one may readily show for a system in full line: 



{lntr{T(x,y,Ai)}, In tr{T(x, y, A 2 )}} = 



(2.6) 



(2.7) 



i.e. the system is integrable, and the charges in involution -local integrals of motion- may 
be obtained by expanding the object In tr{T(x, y, A)}. 

The classical r-matrix associated to the A n affine Toda field theory in particular is given 



n+l 



sinh(A) 



&ii ~~\~ 



i=l 



smh(A) ^ 



(2.8) 



3 Notice that the r-matrix employed here is in fact r^* 2 with r\i being the matrix used e.g. in 
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Note that the classical r-matrix (12.81) is written in the so called principal gradation as is also 
in [IU E] • To express the r-matrix in the homogeneous gradation one implements a simple 
gauge transformation: 

n+l 

r W(A) = V(A) r (rt (A) V(-A) where V(A) = Y j e^^e jj . (2.9) 

Our main aim as mentioned upon is to study the An ATFT on the interval. For this 
purpose we shall employ Sklyanin's formulation (see also [33[ G2] for classical models with 
integrable boundary conditions). It will be convenient for our purposes here to introduce 
some useful notation: 

r ab (X) = r ba (X) forSP, r ab (X) = r*f(A) for SNP 

r* ab (X) = r o6 (A) forSP, r* ab (X) = r£(-A) for SNP 

KtW = n a (X) forSP, r* ab (\) = r t : b (-\) for SNP 

T(A) = T~ 1 (-A) forSP, f(A) = T*(-A) for SNP. (2.10) 

In the situation where non-trivial integrable boundary conditions are implemented one derives 
two types of 'monodromy' matrices T, which respectively represent the classical versions of 
the reflection algebra R, and the twisted Yangian T written in the compact form below (see 
e.g. ©ST]): 

{^(AO, T 2 (A 2 )} = n 2 (Ai - A 2 )T 1 (A 1 )T 2 (A 2 ) - T 1 (A 1 )T 2 (A 2 )f 12 (A 1 - A 2 ) 
+T 1 (A 1 )f* 2 (A 1 + A 2 )T 2 (A 2 ) - T 2 (A 2 )r 1 * 2 (A 1 + A 2 )T 1 (A 1 ). (2.11) 

The modified 'monodromy' matrices, compatible with the corresponding algebras M, T are 
given by the following expressions [2j [11] : 

T(x, y, t, A) = T(x, y, t, A) K~(X) f(x, y, t, A) (2.12) 

and the generating function of the involutive quantities is defined as 

t(x, y, t, A) = tr{K + {\) T(x, y, t, A)} (2.13) 

where c-number representations of the algebra M (T) satisfying ( 12. lip for SP and SNP 
respectively, and also 

{^(AO, Kf{\ 2 )} =0. (2.14) 
Due to (12. lip it can be shown that 

{t(x,y,t,Ai), t(x,y,t,\ 2 )} = 0, Ai, A 2 6 C. (2.15) 
Technical details on the proof of classical integrability are provided e.g. in [2l [Til 134"] . 
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2.1 Classical integrals of motion in ATFT 



We shall exemplify our investigation using the first non-trivial model of the ATFT hierarchy 

DOtl 

for a generic An theory [8]: 



that exhibits both types of boundary conditions, that is the A% case. Recall the Lax pair 



V(x, t, u) = £ d x <S> ■ H + ^ (u e^ H E + e~^ H - - e~^ H E 
2 4V u 

U(x, t,u) = £n-H + j (u e^ H E + e~^ H + ± e~^ H E. e^ H ) (2.16) 

2A 

$, n are n- vector fields, with components (pi, 7Tj, i G {l,...,n}, n = e"+' is the multi- 
plicative spectral parameter. To compare with the notation used in [11] we set ^ = ^- (m 
denotes the mass in [TT]). Note that eventually in [11] both (3, rh are set equal to unit. 

Also define: 

n+l n+1 

E + = ^i? Qi , i?„ = ^i?_ Qi (2.17) 

i=l i=l 

«i are the simple roots, H (n-vector) and i?± Qi are the algebra generators in the Cartan-Weyl 
basis corresponding to simple roots, and they satisfy the Lie algebra relations: 



H, E± a . 



±a>iE± c 



E ai , E.J =^a t -H (2.18) 
J at 

Explicit expressions on the simple roots and the Cartan generators are presented in Appendix 
A. Notice that the Lax pair has the following behavior: 

V*(x,t, -u~ l ) = V(x,t,u), U*(x,t,« -1 ) =U(x,t,u) (2.19) 

where * denotes usual transposition. 

Our objective as mentioned is to examine the system with non-trivial boundaries, thus 
we consider representations of the associated underlying algebras expressed by T. To recover 
the local integrals of motion of the considered system we shall follow the quite standard 
procedure and expand In t(u) in powers of An alternative strategy would be to derive 
the modified Lax pair, compatible with the boundary conditions chosen, and hence the 
associated equations of motion (see e.g. [TT] ) . A systematic derivation of boundary Lax 
pairs independently of the choice of model is discussed in [32] • To expand the open transfer 
matrix and derive the local integrals of motion we shall need the expansions of T(x,y,u), 
T(x,y,u~ l ) and K ± (u). In what follows in the present section we basically introduce the 
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necessary preliminaries for such a derivation, and we also reproduce the known integrals of 
motion for the ATFT on the full line. 

Let T'(x,y,u) = T(x,y,-u~ 1 ) and \J'(x,u) = U(x,-u~ 1 ). Following the logic described in 
[32] for the sine-Gordon model, we aim at expressing the part associated to E + , E_ in U, U' 
respectively independently of the fields, after applying a suitable gauge transformation. More 
precisely, consider the following gauge transformation: 



T(x, y, u) = fi(x) T(x,y,u)Q 1 (y), 

T\x, y, u) = tt-^x) f'(x, y, u) Q(y) tt{x) = e^ x > H . (2.20) 

Then from equation (12.11) we obtain that the gauge transformed operators U, U' can be 
expressed as: 

U(x, t, u) = Q^{x) U(x, t, u) nix) - Q-\x) ^j^- 

dx 

U'(x, t, u) = Q(x) U'{x, t, u) Q-\x) - Q{x) K } . (2.21) 

(JLJb 

After implementing the gauge transformations the operators U, U' take the following simple 
form: 

IJ(x, t,u) = -e-H + —(uE + + -X\ , U'(x, t,u) = ^e-H + ^(uE_ + -X + )(2.22) 
2 4 V u J 2 4 V u / 

where we define: 

9 = n - d x $, 9 = n + d x ®, X_ = e~^ H E_ e^ H , X + = e^ H E + e^*f2.23) 

T, U also satisfy (12.11) . and G, Q are n vectors with components 9i, 9i respectively. 
Consider now the following ansatz for T, T' as \u\ — ► oo [32] 

f{x,y,u) = {I + W{x,u)) exp[Z(x,y,u)] (1+ W{y, u))~\ 

f'(x, y, «) = (!+ W(x, u)) exp[Z(x, y, u)] (I + W(y, u))~\ (2.24) 

where W, W are off diagonal matrices i.e. W = WijEij, and Z, Z are purely diagonal 
z = 1lh=i z uEii- Also 

k=-l k=0 

Inserting the latter expressions ( 12. 25ft in (12. ip one may identify the coefficients and 
Indeed from (12.11) we obtain the following fundamental relations: 

— = v {D) + (u (0) w)^ 

dx 

+ wfj(D) _ fj(D) W + WOJ^W)^ - UC°) - OJ^W)^ = (2.26) 

dx 
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where the superscripts O, D denote off-diagonal and diagonal part respectively. Similar 
relations may be obtained for Z, W, in this case U — > U'. We omit writing these equations 
here for brevity. 

It will be useful in what follows to introduce some notation: 



7^<3>-H = diag(a, b, c), 1 9 • H = diag(a, 6, c), e^"* = 7i 



(2.27) 



explicit expression of a, 6, c and 7, can be found in Appendix B (1B.4[) ; notice that a+6+c = 0. 
From the first of equations (12.261) we may derive the matrices Z, Z . Indeed one may easily 
show that: 



dZ^ 
dx 

dx 



rn 

7 



m 

7 



( -Wff + C'c 



\ 



w£ + (b 



w$ + (b 



\ 



w£ + Ccj 



(2.28) 



it is clear that the latter quantities are zero because of the form of wjp , Wjp see Appendix 



B. Also the higher order Z^ k \ Z^ k ' are given by: 



dx 

dx 
k > 0. 



m 

7 



m 
7 



\ 



(fc-1) 



w& ,/ > + 72 w. 



>,2 



Wt 1] - lsWt 1] J 



(2.29) 



The computation of W, W is essential for defining the diagonal elements. First it is important 
to discuss the leading contribution of the above quantities as \u\ — > 00. To achieve this we 
shall need the explicit form of Z^~ l \ Z^ 1 ^: 



Z^(x,y) 



m{x — y) 
4 



/ vn 
e 3 



e 3 



1 ) 



, Z^\x,y) 



m(x — y) 
4 



/ 



e 3 



\ 



e 3 



1 / 



(2.30) 



The information above will be extensively used in what follows. 
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Before we proceed with the analysis of integrable boundary conditions in ATFT let us 
first reproduce the known local integrals of motion in the periodic case, emerging from the 
expansion — > oo) 

In [trT{u)] = In [tr{(l + W{L, u)) e z(L '~ L ' u] (1 + W(-L, u))" 1 }]. (2.31) 

Notice that in the case of periodic boundary conditions we put our system in the 'whole' 
line (x = L, y = —L), and consider Schwartz boundary conditions, i.e. the fields and their 
derivatives vanish at the end points ±L. Bearing in mind that as u — > — oo the leading 
contribution of e z , (e z ) (see ( 12.301) ) comes from the e Z33 , (e^ 33 ) term, the expression above 
becomes 

7 (k) 

In [ t rT(u - -oo)] = ^ (2.32) 

fc=-i M 

To reproduce the familiar local integrals of motion we shall need both Z(L, —L, u), Z(L, —L, u). 
Let 

I. = ~zg\L, -L, U ) = f <fc( £ <? + £ t , 

P •/ — 1/ i=1 P i=1 



/3 



2 



3l7l f / 8 ?7i 

J 2 = -^Z^{L, —I/, w) = J dx\qp(abc - be) - ^jpbfic + 72« + Isb) 

= 2^3^33 aA-^w) = y da;^— (a&c + fec) - 7^3(7^ + 72« + 7s^ 
(higher local integrals of motion) (2.33) 

the momentum and Hamiltonian (and the higher conserved quantities) of the ATFT are given 
by: 

! rL 2 



I dx^2 (t< $ - Ti 0i) 

i=i 

8=1 ^ 8=1 



2 _A 

/for* 



P 2 = ^(X-2-X 2 ) 

1 f 8 8 * ?7i 2 » \ 

= 2 J L dx {p3^ d ~ abc "> + Tjs ( bc ' + + ^7^(Ti(c - c) + 72(0 - a) + 73(6 - S))J 

^ 2 = ^(J 2 +J_ 2 ) 



1 f / 8 8 * 7T7? * \ 

= g y dx[—(abc + abc) - ~^{bc' - be) - ^(7i(c + c) + 72(0 + a) + 73(6 + S))J 

(2.34) 

Note that the boundary terms are absent in the expressions above, since we considered 
Schwartz type boundary conditions. Also, in the generic situation, for any An \ the sum in 
the momentum V\ and the kinetic term of the Hamiltonian Tii runs from 1 to n, whereas 
the sum in the potential term of the Hamiltonian runs from 1 to n + 1. 



3 SNP boundary conditions 

We turn now to our main concern, which is the study of integrable boundary conditions in 
ATFT. We shall first discuss the boundary conditions that already have been analyzed in 
[TTj . Based on the underlying algebra, that is the classical analogue of the q- twisted Yangian 
we shall reproduce the previously known results [11], so this section serves basically as a 
warm up exercise. In the subsequent section we shall analyze in detail the novel boundary 
conditions (SP) associated to the classical version of the reflection algebra. 

To obtain the relevant local integrals of motion we shall expand the following object 
(consider now x — 0, y — —L): 

In t(u) = In tr^K + (u) T(u) K~(u) T^u" 1 )} 

= In tr^K + (u) fi(0) f{u) Vt~\-L) K~(u) Q{-L) f\u^) fi _1 (0)} (3.1) 

For simplicity here, but without really losing generality we consider Schwartz boundary 
conditions at the boundary point — L and K~(u) oc I. Also K + (u) = i^'(w _1 ) where K is 
any c-number solution of the twisted Yangian. Taking also into account the ansatz ( 12.24ft we 
conclude 

In t(u)=\n trUl + W^O,!!)^- 1 ^) K + (u) fi(0) (1 + W(0,u)) e m-L,n)+&{o,-L,u)\ 

(3.2) 

Recall from the previous section that as u — > — 00 the leading contribution of e z , e z comes 
from the e^ 33 , e Zs3 terms (see (12.301) ). hence 



In t(u) = Z 33 (0, —L, u) + Z 33 (0, —L, u) + ln[(l + W\0, ^fT^O) K + {u) fi(0) (1 + W(0, u)) 

00 7^ _L 7^ 00 f 

E Z 33 + ^33 1 }±_ 

(3 



ir * — ' ir 

k=-l k=0 
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To obtain the explicit form of the boundary contributions to the integrals of motion we 
should first review known results on the solution of the reflection equation for SNP boundary 
conditions. The generic solution for the An case in the principal gradation are given by 
[13122]: 

n+l 

K(X) = (ge x + ge~ x ) E e " + E fv^'^^ + £ ^e^"^"^ 

i=l i>j i<j 

_1 , n+l _ 1_2±I „ _n + l n + l . 

= q 2 4 5 = ±v 2 4 , fij = ±q 4 , = ? 4 , z < j- (3.4) 

In order to effectively compare with the results of [17] as well as being compatible with [11] 
we always express in the text both r and K matrices in the principal gradation. Nevertheless, 
to obtain the matrix in the homogenous gradation as given in [22] we implement the following 
gauge transformation 

K^ h) = V(A) K {p \\) V(-A). (3.5) 

We shall now focus on the A^ case, which is our main example here. Recall that K + (u) = 
-fT^M -1 ) then the iT + -matrix is 3 x 3 matrix written explicitly as: 

K + (u) = u?G + u?F + u~iF + u~?G where 
G = gl, G = gl, 

F = fu + /23 632 + /31 ^13, 

F = / 2 i ei2 + /32 e 23 + /13 e 3 i (3.6) 

and the coefficients g, g, are given in (13. 4p with n = 2. Bearing in mind the explicit 
form of the boundary matrix we may identify the factors fj in the expansion (13. 3p which are 
reported in Appendix C. Taken into account expressions (13. 3p . (IC.1I) and Z^ , given in 
Appendix B we conclude for the first non-trivial boundary integral of motion: 

nf = -^(zfi' + iff+f,) 

2 3 

<**(£« + 0?) + ^ £e^ *) + S(/i2^- $(0) + / M e^-*W - /3xe^- 

i=l ^ i=l ^ 

In general for the A$ ATFT the boundary Hamiltonian with SNP boundary conditions will 
have the following from 

/On 2 n+l n+l 

&<£(*? + ^) + ^ £ e/M ) + E c < ef Q! " (0) ' ( 3 - 8 ) 
L %=1 P i=l i=l 
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which as expected coincides with the boundary Hamiltonian deduced in [llj. It is quite easy 
to check that in the case of a trivial boundary conditions, i.e. K + oc I the boundary terms 
containing c, disappear and the entailed Hamiltonian has exactly the same structure as in 
the bulk case. 

The second conserved charge of the hierarchy is given by 



3m 2 

0,0 o ™2 



1 / / 8 8 771 

2 J L dx \^ abc + - Tjs - ^') - J^^ c + £ ) + 72 ^ a + a ) + 73 ^ + ^ 

2 1 



(/aie-^ + /^e-i" 2 '^ ) - f 13 e~^ ^ + ^(c 2 (0) ~ &(0)c(0) - 6(0)6(0)) 



^ ( /i 2 (c(0) + c(0))ef-- $ (°) - / 23 6(0)ei^-^°) + / a S(0)^ 



9(3 2 
3r? 

^ 2 V 3^ 



— ( " ^(/i 2 e^' $(0) + / 23 e^°> - / 31 e^^(0)) + | (c(0) _ 6(0) + g (0) - &(Q))Y 



(3.9) 

Again when assuming the simplest boundary conditions K + oc I we conclude that all the 
boundary terms containing the factors disappear, exactly as it happens in the first Hamil- 
tonian. The bulk parts of the boundary Hamiltonians above coincide with the ones found 
in the previous section -for Schwartz type boundary conditions. Extra boundary terms are 
added due to the presence of the non-trivial i^-matrix. 

Notice that the boundary analogues of Vk are not conserved quantities anymore similarly 
to the sine-Gordon model on the half line, where only the 'half of the bulk charges are 
conserved after the implementation of consistent integrable boundary conditions. We should 
stress that this is a consequence of the particular choice of boundary conditions, and this will 
become apparent in the next section while analyzing the novel boundary conditions. Note 
also that in the expressions for the boundary Hamiltonian written above there exist no free 
boundary parameter, contrary to the SP case as will see subsequently. Analogous results may 
be seen in the context of quantum integrable spin chains regarding the explicit expression of 
the corresponding Hamiltonians as well as their symmetries [3TL l2T| [22]. 

Let us finally mention that one can in general consider 'dynamical' boundary conditions 
(see e.g. [12], H3j EH]). In this case instead of assuming a c-number solution of the classical 
version of the g-twisted Yangian (12. lip we consider a generic -dynamical- representation of 
the algebra defined as [2]: 

K(A) = L(X - 6) K(X) <g> I L\-\ - 9) (3.10) 

where K is a c-number solution of the classical twisted Yangian [TT1 I22]. and L is any solution 
of the fundamental relation (12.61) e.g. a (/-oscillator. Such boundary conditions for the 
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ATFT have been analyzed in [12]. More precisely, in this case the entries of IK are not c- 
number anymore, but algebraic objects satisfying Poisson commutation relations dictated by 
the underlying classical algebra. At the quantum level these objects, and consequently the 
quantities appearing in the local integrals of motion (13 .7p . (13. 9p . become operators obeying 
commutation relation defined by the g-twisted Yangian. In fact, due to the 'dynamical nature' 
of the boundary conditions extra degrees of freedom, incorporated in L, are attached to the 
boundary. 



4 SP boundary conditions 

We come now to the study of the more intriguing, at least in the present context, boundary 
conditions. Here for the first time we systematically analyze the new boundary conditions 
(SP) starting from the underlying algebra i.e. the reflection algebra. In this case the gener- 
ating function of the local integrals of motion is given by the following expression: 

In t(u) = In tr^K + (u) T(u) K~{u) r -1 (« -1 )} 

= In tr^K + (u) O(0) f(u) Q-\-L) K~{u) Q"\-L) T~ 1 (m~ 1 ) O(0)} (4.1) 

taking into account the ansatz (12.241) we conclude 

In t(u) = In £r{(l + W(0, w)) _1 O(0) K + {u) O(0) (1 + W(0, u)) e z{0 '' L ' u) 

(1 + W(-L, n))- 1 Q- 1 (-L)K-(^)0- 1 (-L)(l + W(-L, w))e-^ (0 '- i ' u) }(4.2) 

The leading contribution of e z , e~ z comes from the e Zl1 , e~ Zl1 terms as iu — > oo, whereas 
as iu — > — oo it comes from the e Z22 , e~ z ' 22 terms. Depending on the limit we assume we 
obtain two distinct expressions for iu — > oo and iu — > — oo respectively: 

In t(iu — > oo) = Zn(0, — L, u) — Zn(0, — L, u) 

+ ln[(l + W(0, it)) _1 O(0) K + {u) O(0) (1 + W(0, u))]u 
+ ln[(l + W(-L, u^Q-'i-L) K~(u) Q'\-L) (1 + W(—L, u))] u 
In t[iu — > -oo) = Z 22 (0, —L, u) - Z 22 (0, -L,u) 

+ ln[(l + W{0, it))- x O(0) K + {u) O(0) (1 + W(0, u))} 22 

+ ln[(l + W(-L, u^Q'^-L) K-(u) 0" X (-L) (1 + W(-L, u))} 22 . 

(4.3) 

Expanding all the terms above we get 

y(k) y{k) OO f+ f 



_k ^ A fc 

k=0 



u k 
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fc=-l fc=0 

Although we follow exactly the same analysis as in the SNP case, we see that the investigation 
of the SP boundary conditions is technically more involved mainly due to the fact that one 
has to consider the behavior of the transfer matrix for both iu — ► oo and —iu — > oo. Another 
technically intriguing point is that the behavior of (1 + W)~ l , which is quite intricate, must 
be studied even if the system is considered on the half line i.e. Schwartz type boundary 
conditions are set at the boundary point — L (see for instance the previous section). 

We shall focus here for simplicity only on diagonal solutions of the reflection equation [26] 
given by the following expressions (in the principal gradation): 

l n 

K (l) (\, £) = sinh(A + iOe- x J2 e ~^ U ~ 1)e n + sinh(-A + i£)e x ^ e^^e^ (4.5) 

3=1 j=l+l 

2 A 

(recall u = e~ ). To obtain the A-matrix in the homogeneous gradation we implement a 
gauge transformation: 

^(A, = V(A) K$(\, f)V(A). (4.6) 

In fact, the presence of non-diagonal boundary conditions does not modify the structure of 
the local integrals of motion, but simply gives rise to more complicated boundary terms. 

Note that in the Ap case we end up with two types of diagonal boundary matrices 
corresponding to the two possible values 1 = 1, 2. We shall consider an example here to 
demonstrate how the particular choice of boundary A-matrix contributes to the integrals of 
motion. Specifically, to obtain the most general results with the least effort it is practical to 
consider a non-trivial left boundary described by A(i), and a right boundary described by 
the A( 2 )-matrix i.e. 

A+K t) = K {1) (u-\ t), K~(u, r) = K&)(% D (4-7) 

The integrals of motion emerging from the first order of the asymptotics of the transfer matrix 
as iu — > ±oo are given by: 

h = z£-z£+hi + hi=-^-(v[ b) - i V3H?) (4.8) 

12m 

(expressions for Z, Z, fj , hf are provided in Appendix B). The momentum and energy are 
directly obtained from the above conserved quantities and defined as: 

-0 2 2 



V[ b) = f dxX;(7r 4 ^-7^^)+X;7r i (O)0 i (O) + |a a .n(O) 
i=i i=i P 



12m -2ij+ -Pa 3 M0) 
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- $>(-L) U-L) - ■ U(-L) + ™le-*-e-*»*W 

1 = 1 

8=1 P 1=1 P " 

Notice the presence of the free boundary parameters in the local integrals of motion above, 
as opposed to the SNP case where no free boundary parameters appear in the corresponding 
integrals of motion. Naturally the two boundary cases are qualitatively distinguished; in SNP 
the c-number iT-matrix contains no free parameters, and consequently no free parameters 
occur in the entailed integrals of motion. In the SP case however the X-matrix contains free 
parameters, which explicitly appear in the boundary integrals of motion. The implementa- 
tion of non-diagonal K-matrices would lead to the appearance of extra boundary terms and 
parameters in the induced local integrals of motion. 

The integrals of motion emerging from the second order of the expansion are derived as: 
h = Z^-Z^+hi + h^^(V?-iV3H?) (4.10) 



where 



Vf = ~ J dx^{abc - abc) + ^{bc' + be 1 ) + ^7i(c - c) + l 2 {a - a) + 73 (& - &)) 
m 2 ( ,„ s ,„ N \ . 3m 2 „„„.a^ 3m 2 



+ 



+ ^e- 2 ^e-^-*(o) ( c(0) + £(0)) + A (&'(0) - b'(0)) + A (& 2 (0) + 6 2 (0)) 
m 2 f , , s , , s \ 3m 2 M 3m 2 



+ ^e- 2 ^e-^(^)( a (_ L)+ a(-L)) + ^(b'(-L)-6'(-L)) +l(6 2 (-L) + 6 2 (-L) 

(4.11 

ft? = ^y"^x(^(a6c + a6c)-| J (6 C / -6cO-^(7i(c + c) + 7 2 (a + a)+73(6 + 6))) 
+ | e -^ e - fe .(o) ( e(0) _ c(0) ) _ !_ ^ (0) + y {0) } + ^ (0) _ g2 (0) 

^c-^e"^-*^) (a(-L) - a(-L)) + -| (tf(-L) + b'(-L)) + 1 (V(-L) - S 2 (-L) 

(4.12) 



4 The parameters £ + , £ are associated to the right left boundary respectively. Note also that there is an 
implicit dependence on the integers . 
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Higher integrals of motions are naturally obtained from the higher order expansion of the 
open transfer matrix but we shall not further pursue this point here. Notice that both Ti.^ 
and Vf are conserved quantities contrary to what happens in the SNP case analyzed in the 
previous section, where only Jx^ are conserved. This is another basic qualitative difference 
between the two types of boundary conditions. Note that from the deduced integrals of 
motion certain sets of equations of motion are entailed. In particular the equations of motion 
arise from the following equations: 

^> = {<>(0,-L), *(«,«)}, ^jM = {«?>«>, -L), *(x,t)}, 
-L<x<0, ie{l,...,n}. (4.13) 

A detailed discussion on the associated equations of motion and the relevant boundary Lax 
pairs systematically constructed along the lines described in [42] will be presented in a forth- 
coming publication. 

As in the analysis of the preceding section for the classical twisted Yangian (SNP) we 
may as well consider dynamical boundary conditions in the SP case. Specifically, one can 
assume a generic -dynamical- representation of the underlying classical reflection algebra 
f l2~H|) defined as 0: 

K(A) = L(X - 6) K(X) ® I L -1 (-A - 0) (4.14) 

where K is a c- number solution of the classical reflection algebra [26], and L is any solution of 
(12. 6j) . Again the extra boundary degrees of freedom are incorporated in L. A more detailed 
analysis of such boundary conditions in the ATFT frame will be presented elsewhere (see 
similar analysis for the sine-Gordon and the vector NLS models in [43] and [34] respectively). 



5 Discussion 

An exhaustive study of the integrable boundary conditions in ArP ATFT was presented by 
systematically deriving the associated local integrals of motion. The key point in our analysis 
is the extraction of the local integrals of motion directly from the transfer matrix asymptotic 
expansion, and there is no conjecture involved as far as their structure is concerned. The 
systematic derivation of the boundary integrals of motion starting from the underlying algebra 
gives rise to two distinct types of boundary conditions associated to the reflection algebra 
and (/-twisted Yangian. 

Noticeably the SP boundary conditions are absent in the analysis presented in [TT] mainly 
because of the a priori strong constraints imposed upon the structure of the boundary con- 
served local quantities. In [11] quantities of the type Vk were a priori disregarded as non 
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conserved -this is true however only for the sine(sinh)-Gordon model (A[ )- whereas as we 
see in the present investigation these objects play a key role in distinguishing the two types 
of boundary conditions! Although sine-Gordon is the prototype model of the class under 
consideration an 'imitation' of its boundary behavior by the higher members of the hierarchy 
could be quite misleading. This is primarily due to the fact that the sine-Gordon is a self- 
conjugate model -soliton and anti-soliton are equivalent entities- and as such it has a very 
peculiar boundary behavior that cannot be naively generalized to higher An ATFT. 

One of the basic differences between the two types of boundary conditions is that in 
the SP case the number of integrals of motion is 'double' compared to the SNP ones. This 
phenomenon not only indicates a qualitatively different behavior of the model as far as the 
boundaries are concerned, but also leads to a modification of the bulk behavior altogether 
(see also e.g. [31]). The 'duplication' of the local integrals of motion in the SP case seems to 
persist to higher orders -we checked explicitly up to third order. More precisely, let Qk, Q-k 
be the local integrals of motion of the An ATFT on the full line, then the boundary conserved 
quantities for each type of boundary conditions are provided by: 

Qt ] = Qk + Q-k + Bk for SNP 

Qf (6) = Q k ±Q- k + B± forSP (5.1) 

Bk, B^ are the relevant boundary terms. In the SNP case only the integrals of motion 
where the bulk part is provided by the sum of Qk, Q-k survive, while in SNP both sums 
and differences provide local conserved quantities, i.e. each one of Q±k (with appropriate 
boundary terms) is a conserved quantity. Moreover in the SNP case no free parameters appear 
in the integrals of motion due to fact that the corresponding c-number K-matrices contain 
no free parameters. However in the SP case, as anticipated, the relevant integrals of motion 
depend on the parameters l ± . It is worth stressing that in the context of integrable 
spin chains the parameters I explicitly appear in the corresponding Hamiltonian as well 
as in the associated symmetry of the model. More precisely, it was shown in [28] that the 
rational open spin chain with diagonal boundary conditions associated to integers is 
gli <S> gln+i-i invariant and U q (gk) <g>U q (gl n +i-i) invariant in the trigonometric case, relevant 
to the ATFT theories. Recall that the lA q {gl n +i) spin chain maybe thought of as an integrable 
lattice version of the An ATFT in the same logic that the critical XXZ spin chain may be 
seen as the lattice version of the sine-Gordon model. 

There exist various studies concerning the underlying symmetry algebras when non-trivial 
integrable boundary conditions are present. Specifically, the symmetry algebra in the context 
of ATFT with SNP boundary conditions -being a twisted algebra- was investigated in [T8] . 
while extensive studies on the underlying algebras in integrable spin chains with both types 
of boundary conditions are presented in [231 121] • An analysis in the spirit of [TBI EH] would 
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provide the non-local integrals of motion forming the exact symmetry algebra in the SP 
case, however this will be presented in a separate publication (see a relevant analysis in the 
quantum case in [23]). 

Another intriguing point associated to the 'folding' of integrals of motion is the possible 
folding of the classical counterparts of Bethe ansatz equations in the SNP case emerging from 
the solution of the spectral problem J32J H5] . Although folding of Bethe anastz equations has 
been reported so far only in isotropic examples we conjecture that it should also occur in 
models associated to trigonometric i?-matrices. In general the structure of Bethe ansatz is 
immediately linked to the underlying algebra, therefore a folding of the associated algebra 
-and the corresponding Dynkin diagrams -would be reflected to the structure of the Bethe 
equations. Extensive studies on the folding of the Bethe equations and the relevant Dynkin 
diagrams are presented in [ZD, [22} [31] . 

In a more physical frame this would be translated to a folding of the associated exact 
boundary S-matrices. Notwithstanding boundary S-matrices were extracted in [28] in the 
SP case, the derivation of boundary S'-matrices in the SNP case is still an open question to 
date in the general case (see e.g. [17]). Having said this the derivation of the Bethe ansatz 
equations for trigonometric spin chains with SNP boundary conditions, and the associated 
boundary S'-matrices will provide significant information at both physical and algebraic level. 

The next natural step would be to identify the relevant boundary Lax pairs for both 
types of boundary conditions along the lines described in [12] . A comparison with the Lax 
pair constructed based on a set of postulates in [11] will be especially illuminating. In the 
SNP case the entailed Lax pair should coincide with that found in [TT] , whereas the Lax pair 
in the SP case will be of a novel from. Generalization of our results for any A n ^ (n > 1) 
ATFT will be also presented in a separate publication. Finally, a similar exhaustive analysis 
regarding principal chiral models (partial results maybe found in [16]) will be particularly 
relevant especially bearing in mind the physical significance of a specific super-symmetric 
principal chiral model within the AdS/CFT correspondence pT7l PfB"] . 

Acknowledgments: I am indebted to J. Avan for useful comments. I wish to thank INFN, 
Bologna Section, and University of Bologna for kind hospitality. This work was partly sup- 
ported by INFN, Bologna section, through grant T012. 

A Appendix 

In this appendix we provide explicit expressions of the simple roots and the Cartan generators 
for An [IS]. The vectors o:, = (a\ , . . . , a") are the simple roots of the Lie algebra of rank 
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n normalized to unity «j • a,, = 1, i.e. 



1 



i-l i+1 



ai— ( , . . . , , 

V ' ' ' V 2i ' V 2i 

Also define the fundamental weights ^ = (/4 



0,...,0 , ie{l,...n} 



(A.1) 



> ■ • ■ j rk j ) 



A*fc) > fc = 1 , ■ ■ ■ , n as (see, e.g., @2]). 



The extended (affine) root a n+ i is provided by the relation 

n+l 

i=i 

We give below the Cartan-Weyl generators in the defining representation: 



(A.2) 



(A.3) 



Eat 
Hi 



i+l 



e»+i 



E n 



i 



n+l 



/'}(<;.> ~ e i+i i+i) ' * = 1 , • • • , n 



(A.4) 



For in particular we have: 

ax = (1, 0), a 2 = 
define also the following 3x3 generators 



1 \/3^ 
2' 2 j ' 



«3 



1 

'2' 2 y 



-El = Eti = ei 2 , i?2 = ^2 = e 23, -E3 = £-3 = —^31 



(A.5) 



(A.6) 



where we define the matrices as (e^-)^ = 5^ 5^. The diagonal Cartan generators H 12 are 

(A.7) 



#1 = ^(en - e 22 ), #2 = ^7|( ei1 + e22 ~~ 2e 3a) 



B Appendix 

From the formulas fl2T26l) . fl2T29|) the matrices Z^, may be determined. In 

particular, we write below explicit expressions of these matrices for the first orders. 

/ ef 1 \ 



w {o) = ^(0) 



eir -1 
yes e 3 o y 



I>(i) 
4 



/ e 2 ^ a c ^ 

-a b 
\ ef c -b / 



o -6 -a \ 

-e~^~b — c 
a — e"3"c y 



(b.i; 
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The higher order quantities are more complicated and we give the matrix entries below for 
W^ 2 \ (define also, C = £): 

= ^(-2 73 + 7i + 72) + j (2a' + b') + ^(-2a 2 - 6c), 

ill >>2 — — /-2 — 

w£> = ^(-273 + 7i + 72) + - c') + ^(c 2 - 06) 

= ^(-2 72 + 7i + 73) + + C) + y (& 2 - oc), 

= ^(2 72 - 71 - 73) + C j(~a> - 2c') + ^(2c 2 + ab), 

Wg ] = -^(2 7 i - 72 - 73) + j (2b' + J) + ^(-26 2 - ac) 

2 ^ 7r >>2 ^ 7r 



3 



^2 j = -^(271 - 72 - 73) + + &') + ^y-(« 2 - 6c) (B.2) 



and 



^i ( 2 2) = ^("272 + 7i + 73) + Q -^(b' - c>) + ^(c 2 - ab), 
wg> = l(-2 l2 + 71 + 73) + y (26' + a') + ^(-26 2 - ac) 
= -\(-2li + 73 + 72) - j (2a' + c') + ^(2a 2 + 6c), 
^3 (2) = ^(27i - 72 - 73) + ^-{V ~ a') + ^(-6 2 + ac), 
W$ = -^(-273 + 72 + 7i) + J@- *) + J (-« 2 + fe) 
Wff = ^(-273 + 7i + 72) + j(b' + 2c') + ^(-2c 2 - aft) (B.3) 

where the prime denotes derivative with respect to x, also a, 6, c, and 7, are defined in (I2.27P 
and have the following explicit forms: 



P ,Qi , 02 \ , ft 1 ®\ t ®2 s P #2 

a = — ;=), o = — =1, c = =, 

2 l 2 2^ 2 l 2 2v/3 2^3' 

7l = e^, 72 = e /'(-^i+^) j 73 = e ^(-^i-^). (B.4) 
Moreover using the expressions above and (I2.29P we have: 

7/7(1) —iS. /■ _*2I > _±?L 

rfzV e 3 m. , (e 3 / / /\ C e 3 / 2 ,2 2n 

-^T = -3— 4 (7i + 72 + 73) + (a' - J) + ^— (a 2 + 6 2 + c 2 ) 

7/7(1) HE. > HE. > HE. 

dZ no earn. C e3 /,/ /\ C e 3 / 2 ,2 2\ 

= -j^ + 7, + 73) + — (*>' - «-) + — (« 2 + 6 2 + «■) 

<._^ (7l + 72 + 73) _^_,)_^ + 62 + c2) 
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dZ\{ earn (ea ? (e> 2 ~ 2 ~, 

"3T = TI (l1 + 72 + 7s) - T " - ° ' + — [a + o + c) 

1/7(1) -is /• -i^ /■ -ii 

dZ 22 e 3 m. . Qe 3 * (e 3 ~ 

= — T < 71 + 72 + 73 > + — <* ~ c > + — < fl + " + c > 

~(7i + 72 + 73) + - f) - |(& 2 + i 2 + 8 2 ) (B.5) 



dx 3 4 3 6 

Finally we report Z 4 - 2 \ Zj? . 



3 

dZ^ e~"3 
3 

(2) 



(72 - 73 - C V - c 2 ') + C 2 ca' + (71c + 72a + 73b) - ( 2 abc) 
( - 7i + 7 3 - C V - a 2 ') + C 2 afr' + (7ic + 7 2 a + 7s&) - C 2 a6c 



dZ 33 



dx 3 
c^ 2) e-i 



I ( - 7i + 7 2 + - fc 2 ') - C V - (71c + 72a + l3 b) + Cabc) 
dx 3 v - 7i + 7 2 - C 2 (&" - & 2 ') + C 2 ^' + (71c + 72a + l3 b) - eabt) 



d7 {2) pit / 

^ = ^3- ( - y 2 + 73 - C 2 (c ' - 2 2 ') + C 2 ch' + (71a + 72 s + 73&) - C 2 a&c) 

f- = ^ ( " Vi + 7 3 + CV - a 2 ') - C 2 ac' - (71c + 726 + 73&) + C 2 Sfcc) . (B.6) 



dzil ] 



dx 

C Appendix 

We present here the boundary contributions in the expansion of the classical open transfer 
matrix for both types of boundary conditions: 

SNP boundary conditions: Recall that in this case the expansion of the generating func- 
tion of the local integrals of motion is given in (13.31) . After some tedious algebra we obtain 
for the boundary terms: 

f = H3g), h = 3^(e fa3 $(0) /31 - e^ {0) fa - e^ )/ 12 ) + ^ (c(0) - 6(0) + 5(0) - a(0)) 
f 2 = -l(/ aie -f«-*<°> + /^ 

- ^(fi2e% ai m (c(0) + 6(0)) - fae^%{0) + f 31 e^^a(0)) - ^(a(O)c(O) + 6(0)c(0) 

+ 1 ((2 7l (o) - 72 (o) - 73(0)) - mo) + fe'(o)) + c 2 (a 2 (o) + 6 2 (o))^ 

- ~ (^(e^- $(0) /3i - e^^fa - e^^°)/ 12 ) + |(c(0) - 6(0) + £(0) - a(0))) 2 . (C.l) 
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SP boundary conditions: We shall need for our purposes here the asymptotics of K ± as 
\u\ — > oo: 



e 



-2ii+ l 



K + {\u\ -> oo, £+) ~ e 33 -- e n + ^e 22 + 0(u~ 3 ) 

AT-(H^cx), r) ~ e n e 33 + ^e 22 + G(^ 3 ). (C.2) 

Then from the expansion of the boundary terms in (14. 3p . (I4.4p we obtain the following explicit 
quantities: 

f + = h+ = ln[^], ^ = -Ce*S(0) + Ce-Tc(0)-e- a ^e-^) J 

f+ = {^ 2 (0Rt(0) - ^-^i(0)^3 4 (0) + ^i(0)fi33 2 (0)(c(0) + c(0)) 

1/ , ^ C 2 /,/, x ,2/ c 2 (0) a 2 (0) 6 2 (0) 6 2 (0) 

- g (2 72 (o) - 7l (o) - 7s(0)) - L(tf (0 ) - c '(o)) + c 2 ( - -U + + -y- + -U 

+ ^1^— fi n (0)fi33 2 (0)(c(0)- C (0)) --(2 72 (0)- 7 i(0)-73(0)) - |- (fe'(0) - c'(0)) 

, 2 / c 2 (0) a 2 (0) & 2 (0) fe 2 (0)^ 
^ V 6 12 12 4/1 

h+ = -Cef6(0)-e- 2i « + e- /3Q3 ^°) 

h+ = {n^(0)!i 3 t(0) - l^-fi^O^O) + !iL (0)^(0) (c(0) + 6(0)) + 3(71(0) - 72 (0)) 

+ £( & (o)-?(o)-6r(o ) + a'(o)) + c a ^ S(0) ■ ' m - m ■ 62(0) a2(0) c2(0) 



+ i 



6 V w v y y ' y 'J ' V 12 12 12 6 12 12 

Ce" 2 ^ 2 ^ , 1 

6 



^{^— fi 2 i(0)fi 3 - 3 2 (0)(c(0) - 6(0)) + - (273(0) - 7i(0) - 72(0)) 
+ £(S<<0) - ,(0) + m - ,(0)) + + M + ffil _ W + m + M) } ,C3) 

Similar expressions are obtained for f~, h~: 

f - = h -=ln[^L^)], ff = Ce-f^-L) _ e -*r e -/^H0 

i- = {tf n (-L)n£(-L) - e ——n{ 1 {-L)n^{-L) + ^-_fi 2 1 (-L)fi 33 1 (-L)(a(-L) + a(-L)) 

1 / j _ * , r , , r ,\ , 2 ( a'(-L) 4c'(-L) 2b'(-L)\ 
+ g - Ti(-L) - 7a(-L)) + C 2 ( - - V ^ + T) 

+ c2 / 2c 2 (-L) a 2 (-L) | 2fr 2 (-L) 



+ z 



12 12 12 



V3{ii-_n 2 1 (-L)n 33 2 (-L)(a(-L) - a(-L)) + - (2 72 (-L) - 7 i(-L) - l3 (-L) 
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/ a'(-L) Ac'(-L) 2V(-L) \ 2 ( 2c\-L) a\-L) 2b\-L) ^ 
^ V 18 18 18 / 4 12 12 12 / J 

h- = -Ce"f b{-L) - e -*%- e -0«*M-L) 

h 2 - = [nU-m^i-L) - e ——n{ 1 {-L)n^\-L) + ^-_^ i( _ L)fi - 3 i ( _ L) ^ ( _ L) + a( _ L) ) 

- o 7l(-^J - 72(-£) + C 7 + a a + 



2V v ' v V " V 6 6 6 6 

„/ c 2 (-L) a 2 (-L) fr 2 (-L) c 2 (-L) a 2 (-L) fr 2 (-L) 
^ V 12 12 6 12 12 12 



+ ^V3{^— n 2 n (-L)Q^(-L) ( - a(-L) + a(-L)) + - ( - 2 l3 (-L) + 7 i(-L) + l2 {-L)) 

/ o'(-L) g(-L) a'(-L) y(-LK 
^ V 6 6 6 6 / 



_ c\-L) _ a 2 (-L) & 2 (-L) _ a 2 (-L) _ 6 2 (-L) _ c 2 (— L) \ , 
12 12 6 12 12 12 
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